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Convergence of approximate deconvolution 
models to the filtered Navier-Stokes Equations 

Luigi C. Berselli * Roger Lewandowski''' 



Abstract 



We consider a 3D Approximate Deconvolution Model (ADM) which belongs to 
the class of Large Eddy Simulation (LES) models. We work with periodic boundary 
conditions and the filter that is used to average the fluid equations is the Helmholtz 
Oh I one. We prove existence and uniqueness of what we call a "regular weak" solution 

■ (watj^at) to the model, for any fixed order N G N of deconvolution. Then, we prove 

I that the sequence {(wat, (7Ar)}jvGN converges -in some sense- to a solution of the filtered 

Navier-Stokes equations, as N goes to infinity. This rigorously shows that the class 
of ADM models we consider have the most meaningful approximation property for 
averages of solutions of the Navier-Stokes equations. 
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'nI" : 1 Introduction 

(N 

It is well-known that the Kolmogorov theory predicts that simulating turbulent flows by 
using the Navier-Stokes Equations 



dtu + V ■ (u^u) - uAu + Vp = f, 
X' (1.1) V-u = 0, 

; u(0,x) = uo(x). 

requires TV = 0{Re^/^) degrees of freedom, where Re = ULu~^ denotes the Reynolds 
number and U and L are a typical velocity and length, respectively. This number TV is 
too large, in comparison with memory capacities of actual computers, to perform a Direct 
Numerical Simulation (DNS). Indeed, for realistic flows, such as for instance geophysical 
flows, the Reynolds number is order 10^, yielding TV of order 10^*^.... This is why one aims 
at computing at least the "mean values" of the flows fields, which are the velocity field 
u = (it^, it^, it^) and the scalar pressure field p. This is heuristically motivated from the 
fact that some gross characteristics of the flow behave in a more orderly manner. In the 
spirit of the work started probably with Reynolds, this correspond in finding a suitable 
computational decomposition 

u = u -|- u' and p = p + p', 



*Dipartimento di Matematica Applicata Universita di Pisa, Via F. Buonarroti 1/c, 1-56127, ITALY, 
berselli@dma.unipi.it, http://www2.ing.unipi.it/~d9378. 

tlRMAR, UMR 6625, Universite Rennes 1, Campus Beaulieu, 35042 Rennes cedex FRANCE; 
Roger.Lewandowski@univ-rennesl.fr, http:/ /perso.univ-rennesl.fr/roger.lewandowski/ 



where the primed variables are fluctuations around the over- lined mean fields. Fluctuations 
can be disregarded since generally in applications knowledge of the mean flow is enough 
to extract relevant information on the fluid motion. 

The "mean values" can be defined in several ways (time or space average, statistical 
averages...); in particular, if one denotes the means fields by u and p, and by assuming 
that the averaging operation commutes with differential operators, one gets the filtered 
Navier-Stokes equations 

dtU + V ■ (u® u) - i/Au + = f, 

(1.2) V-u = 0, 

u(0,x) = u^(x). 

This raises the question of the interior closure problem, that is the modeling of the ten- 
sor R{u) = u (^i u in terms of the filtered variables (u,p). Classical Large Eddy Sim- 
ulations (LES) models approximate i? by w ® w — z^7^(k/kc)D(w) where w ~ u, and 
p(w) = (l/2)(Vw + Vw^). Here is an eddy viscosity based on a "cut frequency" kc 
(see a general setting in [21]). We introduce the new variable w since when using any 
approximation for i?(u), one does not write the differential equation satisfied by u, but an 
equation satisfied by another field w which is hopefully close enough to u. 
Another way, that avoids eddy viscosities, consists in approaching Rhy a quadratic term of 
the form B{v/, w). J. Leray [16] already used in the 1930s the approximation (with our LES 
notation) i?(w,w) = w^w to get smooth approximations to the Navier-Stokes Equations. 
This approximations yield the recent Leray-alpha fashion models, considered to be LES 
models, and a broad class of close models (see e.g. [4, 8, 2, 10, 17]). In [12, 13], we also 
have studied the model defined by B{w,w) = w (g) w, which has also strict connections 
with scale-similarity models. 

The model we study in this paper, is the Approximate Deconvolution Model (ADM), first 
introduced by Adams and Stolz [22, 1], so far as we know. This model is defined by 
i?(w, w) = Z)jv(w) (g)L'jv(w). Here the operator G is defined thanks to the Helmholtz 
filter (cf. (2.2)-(2.4) below) by G(v) = v, where in the paper G = {I — a^A)"^, and the 
operator Dn has the form Dn = '^n=oO' ~ G*)". Therefore, the initial value problem we 
consider is: 

dtw + V • {DN{y^) (E) -DAr(w)) - i^Aw + = f, 

(1.3) V-w = 0, 

w(0,x) = u^(x), 

and we are working with periodic boundary conditions. We already observed that the 
equations (1.3) are not the equations (1.2) satisfied by u, but we are aimed at consid- 
ering (1.3) as an approximation of (1.2). This is mathematically sound since formally 
Dn —^ a := 1 — a^A, in the limit N +co, hence again formally (1.3) will become 
the filtered Navier-Stokes (1.4). What is more challenging is to understand whether this 
property is true or not, in the sense that one would like to show that as the approximation 
parameter grows, then 

w ^ u and Q ^P- 

One would like to prove that solutions of the model converge to averages of the true quan- 
tities, since, we recall that the main goal of LES as a computational tool is to approximate 
the averages of the fiow, which are the only interesting and computable quantities. To 
this end we want to point out that, beside the technical mathematical difficulties, proving 



results of approximation of single trajectories 



w — > u and Q ^ P, 

is not really in the "rules of the game," because the consistency of the model towards single 
(generally strong or not computable) solutions to the Navier-Stokes equations is not the 
most interesting point. Anyway, such convergence is generally known only for a few models, 
as for instance for many of the alpha- models, as a — > O"*". To support again this point 
of view, observe that generally a > is related to the smallest resolved/resolvable scale, 
hence a = 0{h), where h is the mesh size. Reducing a will mean resolving completely the 
flow, hence performing a DNS instead of a LES. 

To our knowledge such a "we// posedness," i.e. proving that w ^ u, is not known for any 
LES model: To our knowledge there are no results showing (not only formally, but also 
rigorously) that the solution of a LES model (w, q) is close or converges in some sense to 
the averages (u, p). 

To continue the introduction to our new results, we recall that model (1.3) has already 
been considered in [14] where we studied the residual stress. It has also been studied 
in Dunca and Epshteyn [7], where it has been proved the existence of a unique "smooth 
enough" solution for periodic boundary conditions. In [7] it is also shown that, that the 
sequence of models (1.3) goes -in a certain meaning- to the Navier-Stokes Equations when 
a 0^, for N fixed. Notice that the model intensively investigated in [12, 13, 3], where 
B(w,w) = w (^1 w, is the special case = and it is also called simplified-Bardina, since 
it resembles some features of the scale similarity models. 

The main topic of this paper is then to study what happens when N goes to infinity 
in (1.3). We prove that the sequence of models (1.3) converges, in some sense, to the 
averaged Navier-Stokes equations (1.2), when the typical scale of filtration (called a > 0) 
remains fixed and the boundary conditions are the periodic ones. Before analyzing such 
convergence we need to prove existence of smooth enough solutions. To this end we needed 
to completely revisit the approach in [7]. To be more precise, let be the 3D torus and 
let (wjv, gAr), with 

wjvGL2([0,T];//2(T3)3)nL°°([0,T];Fi(T3)3), 
qNeL'{[0,T];W''\T^)nL'/-'{[0,T];W^''/-'{Ts)), 

denote the solution of the ADM model (1.3), where T is a fixed time, that can eventu- 
ally be taken to be equal to oo, assuming f G L'^{[0,T]; (H'^iTsf)') and uq G ^^(Ts)^, 
an assumption that we do throughout the paper. We are able to prove existence (cf. 
Theorem 3.1) in such class and our main result is the following. 

Theorem 1.1. From the sequence {(wtv, 9Ar)}AfeN one can extract a sub-sequence (still 
denoted {{'WN,qN)}NeN) such that 

weakly in L'^{[0,T]; H'^{Tsf) n L'^{[0,T]; {T^f), 
strongly in LP{[0,T]; H'^iTsf), VI < p < +oo, 
weakly in L^{[0,T];W^'\T3) D L^/^{[0,T];W^^^/^{T3)), 

and such that the system 

dtw + V • (AwWAw) - uAw + Vg = f, 
(1.4) V-w = 0, 

w(0,x) = u^(x). 



holds in the distributional sense, where we recall that A = G ^ = I — a^A. Moreover, the 
following energy inequality holds: 

(1.5) --^Pwf + z.||V^wf < (f,^w). 

As a consequence of Theorem 1.1, we deduce that the field (u,p) = (Aw, Aq) is a dissipative 
(Leray-Hopf) solution to the Navier-Stokes Equations (1.1). 

Remark 1.1. // one rewrites system (1.2) in terms of the variables w = u and u = A\i = 

Aw, one obtains exactly the system (1.4). This is not a LES model, since it is just a 
change of variables. The LES modeling comes into the equations with the approximation 
of the operator A by means of the family {Z^atItvgn- 

Plan of the paper. Since the paper deals mainly with the mathematical properties of 
the model, we start in Section 2 by giving a precise definition of our filter through the 
Helmholtz equation and we sketch a reminder of the basic properties of the deconvolution 
operator D]\j. The precise knowledge of the filter is one of the critical points in the analysis 
we will perform. We also claim that, beside some basic knowledge of functional analysis, 
we have been able to simplify the proof in order to employ just the classical energy and 
compactness methods. Roughly speaking, we needed to find the correct multipliers and -at 
least in principle- the proof of the main result should be readable also from practitioners. 
Then, we show in Section 3 an existence and uniqueness result for system (1.3). Even if 
this result has already been obtained by Dunca and Epshteyn [7] , our proof is shorter and 
uses different arguments, useful for proving our main convergence result. Indeed, Dunca 
and Epshteyn proved initially a smart but very technical formula about D]\f m. terms of 
series of (— A)'^, but they did not get uniform estimates in N . This is why their proof 
cannot help for passing to the limit when N goes to infinity. Our first main observation 
in this paper is that one can get very easily an estimate uniform in N for 74^/^d]^^(w), 
that also yields estimates for Dn{w) and w itself, always uniform in A^. The leitmotiv of 
the paper is to prove estimates independent of N. 

Finally, Section 4 is devoted to the proof of Theorem 1.1. To prove it, we use the estimates, 
uniform in N, obtained in the proof of the existence result. Another main ingredient of 
the proof, is the derivation of an estimate for dtDN{wi^). Good estimates for this term 
yield a compactness property (a la Aubin-Lions) for Z)jv(wjv), which allows us to pass 
to the limit in the non linear term. We also note that in our argument we keep control 
of the pressure, since it is needed in some arguments and we do not simply neglect it by 
projecting the equations over divergence-free vector fields. 
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2 General Background 
2.1 Orientation 

This section is devoted first to the definition of the function spaces that we use, next to 
the definition of the filter through the Helmholtz equation, and finally to what we call the 
"deconvolution operator." There is nothing new here that is not already proved in former 



papers. This is why we restrict our-selves to what we need for our display and we skip out 
proofs and technical details. Those details can be proved by standard analysis and the 
reader can check them in several references already quoted in the introduction and also 
quoted below in the text. 



2.2 Function spaces 

In the sequel we will use the customary Lebesgue and Sobolev VF^'^ and VF*'^ = 
spaces. Since we work with periodic boundary conditions we can better characterize the 
divergence-free spaces we need. In fact, the spaces we consider are well-defined by using 
Fourier Series on the 3D torus T3 defined just below. Let be given L G R* = {x G M : 
X > 0}, and define Q := [0, L]^ C M"^. We denote by (61,62,63) the orthonormal basis of 
R^, and by x = (xi,X2,X3) G R^ the standard point in R^. We put T3 := 2tx1? jL. Let 
T3 be the torus defined by T3 = (R^/7^). We use || • || to denote the L^^z) norm and 

associated operator norms. We always impose the zero mean condition / ^dx = on 
<f) = w,p, f, or wq. We define, for a general exponent s > 0, 

H, = |w:T3^R^ wGi?"(T3)^ V • w = 0, J wdx = o|, 

where H%T'^)'' = [H'{Ts)]'', for ah A; G N (If < s < 1 the condition V • w = must be 
understood in a weak sense). 

For w G H3, we can expand the velocity field in a Fourier series 

w(x) = ^ Wke"*"***^'^, where k G Tg*" is the wave-number, 

and the Fourier coefficients are given by 

1 

The magnitude of k is defined by 

k :=\k\ = {\ki\^ + \k2\^ + \k3\}-^ . 

We define the norms by 

where of course ||w||q = ||w|p. The inner products associated to these norms are 

(2.1) (w,v)h, = J2 Ikp'wk-^,:, 

ker* 

where here, and without risk of confusion with the filter defined later, Vk denotes the 
complex conjugate of Vk- This means that if Vk = (^k' '^k' '^k)' then Vk = ('^k'^k'^k)- 
Since we are looking for real valued vectors fields, we have the natural relation, for any 
field denoted by w: 



Wk = t;^ / w(x)e dx. 



Wk = w_k, VkGTo 



Therefore, our space is a closed subset of the space Mg of complex valued functions 



equipped with the Hilbert structure given by (2.1). It can be shown (see e.g. [6]) that when 
s is an integer, ||w||^ = HVwlp. One also can prove that for general s € M, (He,)' = EI_s 
(see in [19]). 



2.3 About the Filter 

We now recall the main properties of the Helmholtz filter. In the sequel a > 0, denotes a 
given fixed number and w G H^. We consider the Stokes-like problem for s > —1: 

— a^Aw + w + Vvr = w in T3, 

2.2) 

V-w = inTs, 
and in addition, / vr dx = to have a uniquely defined pressure. 

It is clear that this problem has a unique solution (w, vr) G tls+2 x H^'^^{T3), for any 
w G H^. We put G(w) = w, ^4 = G~^. Notice that even if we work with real valued 
fields, G = maps more generally Mg onto ]H[s_|_2. Observe also that -in terms of Fourier 
series- when one inserts in (2.2) 

w = ^ Wk e^*^ '', 

one easily gets, by searching (w, vr) in terms of Fourier Series, that 

(2.3) w(x)=y- L^Wke+^'^-^ = G(w), and vr = 0. 

With a slight abuse of notation, for a scalar function x we still denote by x the solution 
of the pure Helmholtz problem 

(2.4) Ax = -a^Ax + X = X m T3, G{x) =X- 

and of course there are not vanishing-mean conditions to be imposed for such cases. This 
notation -which is nevertheless historical- is motivated from the fact that in the peri- 
odic setting and for divergence-free vector fields the Stokes filter (2.2) is exactly the same 
as (2.4). Observe in particular that in the LES model (1.3) and in the filtered equa- 
tions (1.2)-(1.4), the symbol " " denotes the pure Helmholtz filter, applied component- 
by-component to the tensor fields Dj\f{w) ^ Di\f{w), u (8) u, and Aw respectively. 



2.4 The deconvolution operator 

We start this section with a useful definition, which we shall use several times in the sequel 
to understand the relevant properties of the LES model. 

Definition 2.1. Let K be an operator acting on H^. Assume that e~*k-x ^j,g eigen-vectors 
of K with corresponding eigenvalues i^k- Then we shall say that ETk is the symbol of K . 



The deconvolution operator is constructed thanks to the Van-Cittert algorithm, and 
is formally defined by 



N 

(2.5) D^:=Y^{l-Gr. 

n=0 

The reader will find a complete description and analysis of the Van-Cittert Algorithm and 
its variants in [f8]. Here we just report the properties we only need for the description of 
the model. 

Starting from (2.5), we can express the deconvolution operator in terms of Fourier Series 
by the formula 

DnM = Yl 5,v(k)wke+^''■^ 

where 

^ / 2|i,|2 \ / 2|i,|2 \ N+l 

(2.6) 5^(k) = ^(3^^j =(l + a^|k|>^,., W = l-(^f^j . 

n=0 \ I I / \ I I / 

The symbol Di\i(k.) of the operator satisfies the following 

(2.7) for each k G 7^ fixed -DAr(k) ^ 1 + a^|kp, as ^ +oo, 

even if not uniformly in k. This means that {DjyjAfGN converges to A in some sense (see 
Lemma 2.2 below). We need to specify this convergence in order to pass to the limit more 
than in "a formal way," to go from (1.3) to (1.4). A general goal for the all paper is to 
precisely determine the notion of D]\f ^ A and to obtain enough estimates on the solution 
w of (1.3) in order to perform such limit. 

The basic properties satisfied hy Dj^ that we will need are summarized in the following 
lemma. 

Lemma 2.1. For each G N the operator ■ Hg — > is self- adjoint, it commutes 
with differentiation, and the following properties hold true: 



(2.8) 1 < 5jv(k) < iV+ 1 VkGTg, 

(2.9) 5^(k) « (Ar + l)l±^^ for large \k\, 

(2.10) lim DnO^) = N + 1 for fixed a > 0, 

|kH+oo 

(2.11) 5jv(k) < (1 + a2|k|2) VkGTa, a>0. 



All these claims are "obvious," in the sense that they follow from direct inspection of the 
formula (2.6). Nevertheless, they call for some comments. A first observation is that (2.10) 
is a direct consequence of (2.9). This says that the spaces are preserved by the operator 
Djsi. More precisely, for all s > 0, the map 

•w 1-^ Di^{w), 

is an isomorphism which satisfies 

\\Dn\\h.=0{N + 1). 



Moreover, the term Z)jv(w) ® Dfyf{w) in model (1.3) is better than the convective term 
Aw (8) in the classical filtered Navier-Stokes Equations, making a good hope for the 
model to have what we call a unique "regular weak" solution (see Definition 3.1 in the next 
section) which satisfies an energy equality. This follows because the sequence {DArjTveN is 
made of differential operators of zero-order approximating A, which is of the second-order. 
The solutions of (1.3) are stronger than the usual weak (dissipative) Leray's solution: this 
is the good news. The bad news is that high frequency modes are not under control and 
may blow up when one lets A'^ to go to infinity, making very hard the question of the limit 
behavior of the sequence of models (1.3), when A'^ goes to infinity. 

In the same spirit of limiting behavior of Djsf -as a byproduct of (2.7) and (2.9)- it can be 
shown that the sequence of operators {0^}^,=^^ "weakly" converges (more precisely one 
has point-wise convergence) to the operator A. The following result holds true. 

Lemma 2.2. Let s G M and let w E Hs+2. Then 

lim D/v(w) = Aw in H,. 

The proof of this lemma is very close to the one of Lemma 2.5 in [15]. Therefore, we skip 
the details. 

Remark 2.1. Since is self-adjoint and non-negative it is possible to define the frac- 
tional powers D'^, for a > 0. From the previous result we also obtain directly that if 
w G Hs, then 

Va>0, lim A""L>^(w) = w in H,. 

Remark 2.2. The reader may observe that most of the properties satisfied by D]\f are also 
satisfied by the Yosida approximation 

which is very common in the theory of semi-groups or in the calculus of variations. To 
compare the behavior of the two approximations, we write the explicit expression for the 
symbol of the Yosida approximation, with A = 1/A^.' 



Ii/iv(k) = (l+a2|k|2) 



^ l + a^|k| 



iV + l + a2|k|2 



(1-ha |k| )aN,k, (^N,k 



iV + l + a^k 



One can directly compute that the asymptotics are essentially the same as in Lemma 2.1, 
but the Van Cittert operator Dn converges to A much faster than the Yosida approximation 
Aiijq, as N goes to infinity. 



3 Existence results 

As we pointed out in the introduction, in this paper we consider the initial value problem 
for the LES model (1.3). The aim of this section is to prove the existence of a unique 
solution to the system (1.3) for a given A^. In the whole paper, a > is fixed as we 
already have said, and we assume that the data are such that 

(3.1) uoGHo, f GL2([0,r];H_i), 
which naturally yields 

(3.2) u^GH2, f GL2([0,r];Hi). 



We start by defining the notion of what we call a "regular weak" solution to this system. 



Definition 3.1 ("Regular weak" solution). We say that the couple (w, is a "regular 
weak" solution to system (1-3) if and only if the three following items are satisfied: 

1) Regularity 



(3.3) w E L2([0, r]; H2) n C([0, T]; Hi), 



(3.4) atwGL2([0,r];Ho) 

(3.5) qeL\%T]-H\T^)), 



2) Initial data 

(3.6) lim ||w(t, •) - u^IIhi = 0, 



3) Weak Formulation 



(3.7) Vv G L2([0,T];i7XT3 



(3.8) 



/ / dtw-v-f I Dn{w)<^Dn{w) -.Vv+u f I Vw : Vv 

Jo Jts Jo Jt-a Jo Jts 

f-T f f-T 



+ / V<7 • V = / / f • V. 

JO Jj-i Jo JT3 



Almost all terms in (3.8) are obviously well-defined thanks to (3.3)-(3.4)-(3.5), together 
with (3.7). The convective term however, needs to be checked a little bit more carefully. 
To this end, we first recall that D]\j maps onto itself, and the Sobolev embedding 
implies that w G C7([0,r];Hi) C L°°([0, T]; ^^(Tg)^). 

Then, we stiU have Dn{^) G C([0,T];Hi) C L°°{[{),T]-L^{T^f). In particular, 

i?7v(w) D7v(w) G L°°([0,T];L3(T3)3)^ 
Consequently, we have at least 



D^(w)$5l)iv(w) GL°°([0,r];i/^(T3)=^)2 CL~([0,T] x T3 



which yields the integrability of -DAr(w) ® Dn{'w) : Vv for any v G L^([0, T]; H^{Ts)^). 

Remark 3.1. We point out that we use the name "regular weak" solution, since (w,g) 
is a solution in the sense of distributions, but it will turn out to be smooth enough to be 
uniqueness and to satisfy an energy equality in place of only an energy inequality, such as 
in the usual Navier-Stokes equations. 

Theorem 3.1. Assume that (3.1) holds, a > and N £ N are given and fixed. Then 
Problem (1.3) has a unique regular weak solution. 

Proof. We use the usual Galerkin method (see for instance the basics for the Navier-Stokes 
Equations in [20]). This allows to construct the velocity part of the solution, since the 
equation is projected on a divergence-free vector field space. The pressure is recovered by 
De Rham Theorem at the end of the process, that we divide into five steps: 



Step 1: we start by constructing approximate solutions w^, solving differential 
equations on finite dimensional spaces (see Definition 3.9 below); 

Step 2: we look for bounds on {w^j^gp^ and {dt'Wm}meN, uniform with respect to 
m € N, in suitable spaces. This is obtained by using an energy-type equality satisfied 
by AV2z)V2 

(wm). Most of these bounds will result also uniform in N, where TV G N 
is the index related to the order of deconvolution of the model; 

Step 3: we use the main rules of functional analysis to get compactness properties 
about the sequence {wmjnGN- This will allow us to pass to the limit when m — > oo 
and is fixed, to obtain a solution to the model; 

Step 4: we check the question of the initial data; 

Step 5: we show that the solution we constructed is unique thanks Gronwall's 
lemma. 

Since Step 1 and 3 are very classical, we will only sketch them, as well as Step 4 which 
is very close from what has already been done in [5, 19, 24]. On the other hand. Step 2 
is one of the main original contributions in the paper and will also be useful in the next 
section. Indeed, we obtain many estimates, uniform in N, that allow us in passing to 
the limit when N goes to infinity and proving Theorem 1.1. Also Step 4 needs some 
application of classical tools in a way that is less standard than usual. We also point out 
that Theorem 3.1 greatly improves the corresponding existence result in [7] and it is not 
a simple restatement of those results. 

Step 1 : construction of velocity's approximations . 

Let be given m G N* and define to be the space of real valued trigonometric polynomial 
vector fields of degree less or equal than n, with vanishing both divergence and mean value 
on the torus T3, 

(3.9) := {w G Hi : [ w(x) e'"'^'' = 0, Vk, with |k| > m}. 

The space has finite dimension, denoted by dm- Moreover, C V^+i and, in the 
meaning of Hilbert spaces. Hi = UmgN*^^. We notice that is a subset of the finite 
dimensional space 

:= {w : T3 - C3, w = ^ Wk e+'''-^}, 

keT3,|k|<m 

and the space Vm can be described as 

(3.10) V^:=W„nHo. 

We denote by (ei, . . . , e^^^^) an orthogonal basis of V^- Remark that this basis is not made 
of the e"'"*'''-^'s. Nevertheless, we do not need to know it precisely. Moreover, the family 
{ejljgN is an orthogonal basis of Hq as well as of Hi. As we shall see in the following, 
the ej's can be chosen to be eigen- vectors of A, with ||ej|| = 1. 

Let Pm denote the orthogonal projection from H^ (s = 0, 1) onto V^- For instance, for 
Wo = = I]^i w'jej, we have 

dm, 



In order to use classical tools for systems of ordinary differential equations, we approximate 
the external force by means of a standard Friederichs mollifier, see e.g. [23, 24]. Let p be an 
even function such that p € Co"(IR), < p{s) < 1, p{s) = for \s\ > 1, and f^p{s) ds = 1. 
Then, set F(t) = f(t) if t € [0,T] and zero elsewhere and for all positive e define f^, the 
smooth (with respect to time) approximation of f, by 

lit) :=i /" pf^) F{s)ds. 

Well known results imply that if (3.1) is satisfied, then ^ f in L^([0, T]; Hi). Thanks 
to Cauchy-Lipschitz Theorem, we know the existence of a unique function 

i=i 

and of a positive such that the vector {wm,i{t), ■ ■ ■ ,Wm,drn{'t)) is a solution on 
[0,Tm] ^ [0,T], with WmjiO) — w^ j, in the sense that Vv G V^, Vt € [0, T^] it holds 

(3.11) 

/ i9tWm(i,x) • v(x) dx- {DN{-Wrn) «) DN{'Wm)){t,x) : Vv(x) dx 

-u / Vwm(t,x) : Vv(x) dx = / f i/„(t, x) • v(x) dx, 
7T3 JT'a 



where 

As we shall see it in step 2, we can take = T. This ends the local-in-time construction 
of the approximate solutions Wm(t,x). □ 

Remark 3.2. We want to stress to the reader's attention that a more precise notation 
would be 

instead of Wm . We are asking for this simplification to avoid a too heavy notation, since 
in this section both N and a are fixed. 

Step 2. Estimates . 

As in the classical Galerkin method, we need some a priori estimates, first to show that 
the solution of the {dm x (im)-systems of ordinary differential equations satisfied by Wmj 
exists in some non-vanishing time-interval, not depending on m (to this end a energy-type 
estimate is enough). Next, we want to obtain estimates on the w^'s and the dtWm for 
compactness properties, to pass to the limit, when m — > oo and N is still kept fixed. 
As usual, we need to identify suitable test vector fields in (3.11) such that, the scalar 
product with the nonlinear term vanishes (if such choice does exist). We observe that 
the natural candidate is ADi\f{'Wm)- Indeed, since A is self-adjoint and commutes with 
differential operators, it holds: 



{DNi^m) DN{y^ni)) ■ V{ADN{Wm)) dx 

G[DN{^m) ® DN(Wm)) ■ V(ADAr(Wm)) dx 
[ iAG){DNiWm) ® DNi^m)) : V {D N (Wm)) dx = 0, 



because Ao G = ld on H^, V • [Dj\[{w„i)) = 0, and thanks to the periodicity. This yields 
the equahty 

(3.12) {dtWm,ADNiy^m)) - l^{AWm, ADNiWm)) = (f l/m , ^^Af (Wm)) ■ 

This formal computation asks for two clarifications: 

i) We must check that AL)jv(wm) is a "legal" test function, to justify the above formal 
procedure. This means that for any fixed time t, we must prove that ADN{wm) G 

ii) Estimate (3.12) does not give a direct information about Wm itself and/or dt^m- 
Therefore one must find how to deduce suitable estimates from it. 



Point i) is the most simple to handle. On one hand, we already know that G(Ho) = 
H2 C Hq. On the other hand, formula (2.3) makes sure that G{Wm) C W^. We 
now use representation (3.10), and we deduce that G(Vm) C V^- Finally, it is clear 
that Ker(G) = 0. Therefore, since has a finite dimension, we deduce that G is an 
isomorphism on it. Then the space is stable under the action of the operator A 
as well as under that of D]\j. This makes AD]\]{\^m){t^ ■) & a "legal" multiplier in 
formulation (3.11), for each fixed t. Moreover, since A and Dpf are self-adjoint operator 
that commute, one can choose the basis (ei,-- - ,erf^,---) such that each is still an 
eigen-vector of the operator A and Dj^i together. Therefore, the projection commutes 
with A as well as with all by-products of A^ such as Dn for instance. We shall use this 
remark later in the estimates. 

The next point ii) is not so direct and constitutes the heart of the matter of this paper. 
The key observation is that the following identities hold: 

(3.13) {dt^m,ADN{^^)) = ^^pi^l (w„)||2, 

(3.14) (- Aw„,ADAr(w™)) = \\VA^Dl{^m)f, 

(3.15) (f ADMi^ra)) = {AhDl{ly^),A\Dl{^m)) . 

These equalities are straightforward because A and Dn both commute, as well as they do 
with all differential operators. Therefore (3.12) or better 

i|pll)|(w„)||2 + ||V^^4(w„)||2 = (^5Z)|(f,/„), AiDj,(w„)) 

1 1 

shows that the natural quantity under control is A2Z)^(wm). As we shall see in the 
remainder, norms of this quantity do control w-m, as well as the natural key variable 
-C^7v(wm)- Finally, this yields an estimate for dt^m- 

Since we need to prove many a priori estimates, for the reader's convenience we organize 
the results in the following Table (3.16). We hope that having a bunch of estimates 
collected together will help in understanding the result. In a first reading one can skip the 
proof of the inequalities, in order to get directly into the core of the main result. 
The results are organised as follows: In the first column we have labeled the estimates. 
The second column precises the variable under concern. The third one explains the bound 
in term of space function. The title of the space means that the considered sequence is 
bounded in this space. To be more precise, Em,N € F for any variable Em,N and a space 



F, means that the sequence {Em,N} (jn,N)<m'^ is bounded in the space F. Finally the fourth 
column precises the order in terms of a, m, and N for each bound. Of course each bound 
is of order of magnitude 

C(^||uo||l2 + -||f||L2([o,T];L2)^, 

and this why we do not mention it in the table. We also stress that all bounds are uniform 
in m. All bounds except (3.16-h) are uniform also in A^. Moreover, we shall also see that 
they are uniform in T yielding, Tm = T for each T. We mention that we could take T = oo 
at this level of our analysis, so far the source term f is defined on [0, oo[. 



(3.16) 



Label 


Variable 


bound 


order 


a) 


1 i 


L°°([o,r] 


Ho)nL2([o,r] 


Hi) 


0(1) 


b) 




L°°([0,T] 


Ho)nL2([o,T] 


Hi) 


0(1) 


c) 


Dj!, (wm) 


L°°([0,T] 


Hi)nL2([o,r] 


H2) 


O(a-i) 


d) 




L°°([0,T] 


Ho)nL^([o,r] 


Hi) 


0(1) 


e) 


Wm 


L°°([0,T] 


Hi)nL^([o,r] 


H2) 


O(a-i) 


f) 




L°°([0,T] 


Ho)nL^([o,r] 


Hi) 


0(1) 


9) 




L°°([0,T] 


Hi)nL^([o,r] 


H2) 


0{a-' + {N + 1)) 


h) 


dt^m 


^^([0,T];Ho) 





Checking (3.16-a) — For the simplicity, we shall assume that f G L^{[0,T] x ¥3)^, but 
the proof holds true also for f € -L^([0, T]; H_i): one has to substitute in (3.15) the 

integral over T3 with the duality pairing ( . ) between Hi and H_i and estimate in a 

- 1 1 _ 1 1 

standard way the quantity (^fi/^, ADj\[{wm)) = {A^ D'^{{i/^), D^{-Wm))- When one 

integrates (3.12) with respect to time on the time interval [0,t] for any time t < T^, (by 

using (3.13)-(3.14)-(3.15), and Cauchy-Schwartz inequality) one gets 

(3.17) 

^P^D|(w„)(t,.)f + ||VAiD|(w„)||2dr 

< hA^DlFrnm^ + I plz?J(f 1/^)11 • \\AhDl{^N,r.)\\ dr. 
^ Jo 

1 — — 1 — 1 — 1/2 

Notice that A^D'^{if:) = A^2D^(f^). Since the operator A~'2Dj^ has for symbol p^i^ < 1, 
1 i_ 

then ||j42Z)^f^|| < C||f|| (cf. also Remark 2.1 and the properties of classical mollifiers). 
Since Pm, commutes with A and Djsf we also have 

P5d|^p^u^|| = ||p^^Id|u^|| < wa'^d^u^w < ||uo||. 

By using Poincare's inequality combined with Young's inequality, and standard properties 
of mollifiers, one gets 

(3.18) ipiz)|,(w„)(t, Of + ^ ^* WVA'^DlMf dr < C(||uo||, Mmio^TW.,))- 

When one returns back to the definition of w^, one obtain (as a by product of (3.18) and 

1 1 

also because the e^'s are eigen-vectors for A and D]\j and therefore for A^D^) that 

y^,PN,jWm,j{tf < C(||uo||, ||f||L2([o,T];H_i))i 



making sure that one can take = since no pN,j vanishes. 

Checking (3.16-b)-(3.16-c) — Let v G H2. Then, with obvious notations one has 



\\A-2vf = ^ (1 + a'|kp)|?kP = ||vf + a'||Vv||2. 

1 1 

It suffices to apply this identity to v = D^{wm) and to v = diD^{wm) {i = 1,2,3) 
in (3.17) to get the claimed result. 

Checking (3.16-d)-(3.16-e) — This is a direct consequence of (3.16-b)-(3.16-c) combined 
with (2.8), that we also can understand as 

||w|l, < \\Dn{w)\\s < (iV + l)||w||„ 

for general w and for any s > 0. This explains how important is to have a "lower bound" 
for the operator D^r, since estimates on D^w, a > 0, are in inherited by w. 

Checking (3.16-f) — The operator A^Dj^ has for symbol (1 + a^|kp)p]y^^ while the one of 

1 i 

Dn is (1 + a^|kp)pAr^k- Since < pN,]i < 1, then ||DAf(w)||s < ||^2L)^(vir)||s for general 
w and for any s > 0. Therefore, the estimate (3.16-f) is still a consequence of (3.16-a). 

Checking (3.16-g) — This follows directly from (3.16-e) together with (2.8). This also 
explains why the result depends on N, since we are using now the upper bound on the 
norm of the operator Djsf, while in the previous estimates, we used directly the equation 
as well as the lower bound for the Van Cittert operator D^. 

Remark 3.3. The fact that this estimate is valid for each N , but the bound may grow with 
N is the main source of difficulties in passing to the limit as N ^ +00. Also the lack of this 
uniform bound requires some work to show that (certain sequences of) solutions to (1.3) 
converge to the average of a weak (dissipative) solution of the Navier-Stokes equations. 

Checking (3.16-h) — Let us take 5jWm G as test vector field in (3.11). We get 
where 



(3.19) := V • (L>Ar(w„) (g) DN{Wm)). 

So far WmiO,-) = IPm(uo) £ H2 and obviously ||lPm(uo)||2 < C*" ""^IIuqH, we only have to 
check that Ajy^rn is bounded in L^([0,T] x ¥3)^. Thanks to (3.16-f), it is easy checked 
with classical interpolation inequalities, that Diy{wm) G ^'^([O, T]; L^(T3)^). Therefore, 
-DAr(wm) DN(wm) e L^{[0,T];L^/^{T3)^). Because the operator (V-) o G makes to 
"gain one derivative," we deduce that Atv,^ G L'^{[0,T];W''^^/'^{Tsf), which yields to 
Aiv,m 6 L\[0,T] X ¥3)3 since W^^^/^{Ts) C L^{Ts) C L^T^) and L\[0,T]; L^T^f) is 
isomorphic to L^([0, T] x Ts)^ (see [19]). Moreover, the bound is of order 0{a~^) as well, 
because of the norm of the operator (V-) o G that we do not need to specify, but where 
a is involved. The remainder of the proof is a very classical trick, based on Gronwall's 
lemma (see for instance in [18] for a detailed report about this method). Notice that this 
bound on the nonlinear term is not optimal, but fits with our requirements. □ 



Step 3 : Passing to the limit in the equations when m ^ oo, and N is fixed. 

Thanks to the bounds (3.16) and classical tricks, we can extract from the sequence 
{wmjneN a sub-sequence converging to a w G L°°([0, T]; Hi) n L^([0, T]; H2). Using 
Aubin-Lions Lemma (here one uses (3.16-h) and again classical tricks), one has: 

(3.20) w„ ^ w weakly in L2([0,r];H2), 

(3.21) ^ w strongly in LP ([0,r]; Hi), Vp < 00, 

(3.22) dt^m ^ dtw weakly in ^^([o^ t]; Hq). 

This already implies that w satisfies (3.3)-(3.4). From (3.21) and the continuity of in 
H<j we get that D]\j{wm) converges strongly to Dj^{w) in ^"^([0, T] x T3). Then, L>7v(wm)(g) 
DNi^m) converges strongly to Dj^{w)i^Diy{w) in L^([0, T] x T3). This convergence of w, 
together with the fact that fi/m converges strongly, implies that for all v € L^([0, T];Hi) 



(3.23) 



/ / dtw ■ V d:x. dr — / / L)7v(w) DAf(w) : Vvdxdr 
JO JT3 Jo JT3 ^ ^ 

+v / / Vw : Vvdxdr = / / f-vdxdr. 

Jo JT'j, 



Arguing similarly to [18], we easily get that w satisfies (3.6). 

We are almost in order to introduce the pressure. Before doing this, let us make so 
comments about the variational formulation above. We decided to take test vector fields 
in L^([0, T];Hi) to be in accordance with classical presentations. The regularity of w 
however yields V- {D n (w) D n (w)) G L'^{[0, T] x ¥3)^ as well as Aw € L^{[0, T] x T^f. 
Consequently, one can take vector test fields v € L^([0,T];Ho) in formulation (3.23) that 
we can rephrased as: Vv € L^([0, T]; Hq), 

(3.24) [ [ (9iW + Aat - z^Aw- f) • vdxdr = 0, 

Jo Jts 

where for convenience, we have set 

An :=V ■ (-Div(w) 0L»iv(w)). 

Therefore, for almost every t € [0, T], 

¥{t, •) = {dtw + An- uAu - f)(t, •) e L^Tsf 

is orthogonal to divergence-free vector fields in L^(T3)'^ and De Rham's Theorem applies. 
Before going into technical details, let us first recall that among all available versions of 
this theorem, the most understandable is the one given by L. Tartar in [23], a work that has 
been later reproduced in many other papers. Notice also that there is a very elementary 
proof in the periodic case [19]. Now, from (3.24), we deduce that for each Lebesgue point t 
of F, there is a scalar function q{t, •) G H^{T3), such that F = — Vg, that one can rephrase 
as 

(3.25) dt-w + AN-vAw + Vq = l. 



Without loss of generality, one can assume that V • f = 0. Therefore, taking the divergence 
of equation (3.25) yields 

Aq = V- An, 



which also easily yields q G L^([0, T]; //^(Ts)), closing this part of the construction. □ 



Step 4 : About the initial data. 

We must check that w(0, •) can be defined and that w(0, •) = uq, such as defined in (3.6). 
Thanks to the estimates above (mainly S^w G L^([0, T]; Hq), together with the regularity 
about w that we get) one obviously has w € C([0,T];Hi). This allows us to define 
w(0, •) € Hi and also to guarantee 

lim ||w(t,-)-w(0,-)||Hi =0. 



It remains to identify 'w{Q, •). The construction displayed in Step 1 yields for m € N, 

/o 



(3.26) Wm(t,x) = Pm(uo)(x) + /" 9jWm(s,x)(is, 

JO 



an identity which holds in C^([0, T] x il). Because of the weak convergence of (c?fWm)mgN 
to dt'w in L^([0, T];Ho) and thanks to usual properties of the projection's operator 
one can easily pass to the limit in (3.26) in a weak sense in the space L^([0, T]; Hq), to 
obtain ^ 

w(t, x) = uo(x) + / 9(w(s, x) (is. 

JO 

Therefore, w(0,x) = uo(x), closing the question about the initial data. □ 
Step 5 : Uniqueness. 

Let wi and W2 be two solutions corresponding to the same data (uq, f) and let us define as 
usual W := wi — W2. We want to take ADm{^) as test function in the equation satisfied 
by W, because we suspect that this is the natural multiplier for an energy equality. But 
before doing this, we must first check that this guy lives in L^([0, T] x ¥3)^ to be sure 
that this is a "legal" multiplier. Notice that ADjy has for symbol (1 + a^|kp)^/97v,k ~ 
(iV + l)(l + a2|k|2)2/a2|k|2 ^ (Af + l)a2|k|2 for large |k|. Therefore, for each fixed N G N, 
ADi\f is like a Laplacian and makes us loose two derivatives in space. (Only in the limit 
—>■ +00 we will loose four derivatives!) Fortunately, W G L^([0,T];H2) and therefore 
ADj\f{W) G L'^{[0,T] X Ts)^, making this guy a good candidate to be the multiplier we 
need. After using tricks already introduced in this paper, we get 



\A2D^(W)f + i^\\VA2D^(W) 



1 d 

(3 27) < |((I?7v(W).V)Div(w2),I?^(W))|, 

< ||I)^(W)||2,^^^^||VI)^(W2)||, 

< ||I)^(W)r/2||VZ)^r(W)f/2j|VZ);v(w2)||, 

where the last line is obtained thanks the well-known "Ladyzhenskaya inequality" for inter- 
polation of L'^ with !?■ and H^, see [11, Ch. 1]. Starting from the last line of (3.27), we use 

the inequality ||D7v(W)|| < \\A^d\{W)\\ together with ||D7v(VW)|| < piD|^(VW)||, 
the fact that -Dat and V commute, H-DtvH = (A -|- 1), the bound of W2 in L°°([0, T]; Hi), 
and Young's inequality. We obtain 

i|pi4(w)f + H|VAi4(w)f 



27(A-M)4supi>o ||Vw2|l 1 1, v„ 1 , 

< ^ on.^ - ^P2Z)^(W)||2 + -\\VA-2D%{W) 



In particular, we get 



1 d 1 i o 27(iV + l)'^sup.>n IIVW2II 1 i 



1 



1 1 1 



We deduce from Gronwall's Lemma that A2 D^CW) = because A2 D^{'W){0, ■) = 0. 

1 1 

To conclude that W = 0, we must show that the kernel of the operator A^ is reduced 
to 0. This is trivial, since this operator has for symbol (1 + a^|kp)/jAr^k ~ a|k| for large 
values of k. This symbol never vanishes and the equivalence at infinity shows that this 
operator is of same order of a|V|. Therefore, it is an isomorphism that maps onto 
H^-i. This concludes that the considered kernel is reduced to zero, proving uniqueness of 
the solution. □ 

Remark 3.4. As we have seen, ADjv(w) is a legitimate test function in Equation (3.25). 
When using computation rules already detailed in the paper, we get the following energy 
equality satisfied by A^/'^D^J'^{w), 

l^^\\Ay^DfMf + u\\VAy^DfMf={A~^^^^^^ 

As we shall see in the sequel it seems that it is not possible to pass to the limit N +00 
directly in this "energy equality" and some work to obtain an "energy inequality" is needed. 



4 Passing to the limit when N ^ oo 

The aim of this section is the proof of main result of the paper, namely Theorem 1.1. We 
now denote, for a given G N by (-wj^f^qiy), the "regular weak" solution to Problem 1.3. 
For the sake of completeness and to avoid possible confusion between the Galerkin index 
m and the deconvolution index A^, we write again the system satisfied by wat 



^tWAT + V • (L>Ar(wAr) (g)Z)Ar(wAr)) - i/Awjv + VgAT = f in [0,r]xT3, 
(4.1) V-Wiv = in [0,r] x T3, 

WAr(0,x) = ujj(x) in T3. 

(More precisely, for all G N we set wat = limm^+00 Wm,Ar,a), and the scale a > 
is fixed.) We aim to prove that the sequence {{y^N, Qn)}n€N has a sub-sequence which 
converges to some {w, q) that is a solution to the averaged Navier-Stokes Equations (1.4). 
Recall that {Aw, Aq) will be a distributional solution to the Navier-Stokes Equations. This 
result gives a undeniable theoretical support for the study of this ADM model. 

We divide the proof into two steps: 

1. We search additional estimates, uniform in A^, to get compactness properties about 
the sequences {Dn{wn)}n&i and {wTvjArgN; 

2. We prove strong enough convergence in order to pass to the limit in the equa- 
tion (4.1). 

Of course the challenge in this process is to pass to the limit in the nonlinear term 
I?Ar(w7v) (g) D]y{-W]\f). This is why we seek for an estimate about dtD^iy^N), knowing 
that we already have some estimate for D]\f(\V]\f). This is how we get a compactness 
property satisfied by D]\f{w]\f), that we use for passing to the limit in the nonlinear term. 



Step 1 : Additional estimates . 

We quote in the following table the estimates we shall use for passing to the limit. The 
Table (4.2) is organized as the previous one (3.16). Recall that -for simplicity- G F 
for any variable En and a space F, means that the sequence {E^j^^f^ is bounded in the 
space F. 



Label 


Variable 


bound 


order 


a 


WiV 


L-([0,r];Ho)nL^([0,T];Hi) 


0(1) 


b 


WAT 


L-([0,r];Hi)nL^([0,T];H2) 


O(a-i) 


c 




L-([0,r];Ho)nL^([0,T];Hi) 


0(1) 


d 


dtWN 


L^([0,T] xTa)^ 


0{a-') 


e 


QN 


L^([0,T];i7i(T3))nLV^([0,T];t^^'5/^(T3)) 


O(a-i) 


f 




L4/^([0,r];H_i) 


0(1) 



Estimates (4.2-a), (4.2-b), (4.2-c), and (4.2-d) have already been obtained in the previous 
section. Therefore, we just have to check (4.2-e) and (4.2-f). 

Checking (4.2-e) — As usual, to obtain regularity properties of the pressure (at least in 
the space-periodic case), we take the divergence of (3.25) obtaining 

-AqN = V • Atv - V ■ f , 

where we recall that 

An = V • (DAr(wAr) (g) Dn{wn))- 

Next, since f G L2([0, T] x ¥3)^, then we get V 1 e L'^{[0,T]; H^{Tsf). We now investigate 
the regularity of A^r . We already know from the estimates proved in the previous section 
that Atv € L^([0,r] X Tsf. This yields the first bound in L'^{[0,T]; {T3)) for q^. 

We now seek the other estimate for q^. Classical interpolation inequalities combined 
with (4.2-c) yield Dn{wn) G L10/3([o, T] x T3). Therefore, A^v G L^^^{[0,T];W^'^/^{Ts)). 
Consequently, we obtain 

qN e L\[0,T];HHTs))nL'/^{[0,T];W^''/^{Ts)). 



Checking (4.2-f) ~ Let be given v G ^^([0, T]; Hi). We use D7v(v) G L'^{[0,T];Ui) as 
test function in the equation satisfied by {wiy,qiy), that is now (by the results previously 
proved) a completely justified computation. We get, by using that 9fW G L^([0,T] x T^)'^ 
(as well as all other guys in the equation), commutes with differential operators, G 
and Dn are self-adjoint, the pressure term cancels because V • Diy{v) = 0, and classical 
integrations by parts 

{dtWN,DNiv)) = {dtDN{y^N),v) 

= u{AwN,DN{^r)) - {Dn{wn) ® Z?jv(wjv), A7(W)) - (I?jv(f), v). 

We first observe that 



(4.4) |(Aw^,L>^(v))| = |(VZ)jv(w7v),Vv)| < Ci(t)||v||i, 

and we use the L'^{[0,T]; (T^)^) bound for Diy{-W]\f), to infer that the function Ci{t) G 
L^([0, T]), with a bound uniform with respect to G N. Using ||i57v(f)|| < ||f|| already 



proved in the previous section and Poincare's inequahty, we handle the term concerning 
the external forcing as follows: 

(4.5) |(Z)^(f),v)| <C||f||||v||i, 

C being the Poincare's constant. Finally, from (4.2-c) and usual interpolation inequalities, 
we obtain that Dn{^n) belongs to L^/^{[Q,T]]L'^{T^f), which yields 

Dn{^n) (i?) Dn{^n) L^'\[G,T]-L\T^f). 



Therefore, when we combine the latter estimate with ||L)jv(Vv)|| < ||Vv||, to get 

(4.6) \[Dm{^n) ® Dn{v^n),Dn{^^)) \ < C2(t)i|v||i, 

where C2(t) € L'*/^([0,r]) and it is uniform in G N. The final result is a consequence 
of (4.3) combined with (4.4), (4.5), (4.6), and Ci{t) + ||f(t,-)ll + C2{t) G L^/^{[Q,T]), 
uniformly with respect to G N. □ 

Step 2 : Passing to the limit . 

From the above estimates and classical rules of functional analysis, we can infer that there 
exist 

wGL°°([0,T];Hi)nL2([0,r];H2), 

zGL°°([0,r];Ho)nL2([0,r];Hi), 

gGL2([0,r];i7i(T3))nL5/3([0,T];Ty2,5/3(Tr3)), 

such that, up to sub-sequences, 

r weakly in ^^([O, T]; H2), 

WAT — ^ w I weakly* in L~([0, T]; Hi), 

1^ strongly in LP ([0,r]; Hi), Vp < 00, 

dt^N — > dtw weakly in L^{[0, T] x T3), 

(4.7) r weakly in L2([0,r] ; Hi), 
Dn{wn) — > z I weakly* in L~([0, T]; Hq), 

[ strongly in LP{[0,T] x ¥3)^, Vp < 10/3, 

dtDNi^N) dtz weakly in L4/3([o, T]; H_i), 

QN q weakly in ( [0, T] ; (T3) ) n L^/^ ( [0, T] ; VF^'S/a (T3 ) ) . 

Of course, we have 

(4.8) Dnv^n(^Dnv^n — >z(g>z strongly in LP([0, T] x ¥3)^ Vp < 5/3. 

All other terms in the equation pass easily to the limit as well. Our proof will be complete 
as soon as we shall have checked that z = Aw, thanks to (4.8). However, this is almost 
straightforward as we shall see, the hard job being already done by the proof of the various 
estimates. 

Let us consider v G L^([0, T]; H2). We have (L>7v(wAr), v) = (wat, DAr(v)). We claim that 
-CAf(v) Av strongly in L^([0,T] x ¥3)^. This will suffice to conclude the proof. Indeed, 



if such a convergence result holds, we have using (4.7) and computational rules already 
quoted, 

(i:>Ar(wAr), v) = (w^iv, £'Ar(v)) 



Z, V 



(w, Av) 



Z, V 



(^W, v) 



yielding z = ^w. 

It remains to check the claim. We can write 



V = J] Vk(t) 



ikx 



ker* 



and consequently (see in [19]) 



Il2([o,t] 



H2) = E 1^1' / l^k(t)l'dt < OO. 



Let e > being given. Then, there exists < K = K{v) G N such that 



Since < (1 — j^) < 1, we have 

fT 



[ \\{A- DmM^ = V (1 + a'|k|2)2(l - p^^^f f \^^{t)\^dt, 
Jo kg^^* Jo 

0<|k|<X 
\k\>K 

< Yl (l + «'|kp)'(l-Piv,k)'^'^|vk(t)pdt+ |. 



o<|k|<i<: 

Observe that - for each given k G - we have pn,Il 1 when N oo. Therefore, there 
exists Nq (obviously depending on v and on K) such that for all N > Nq, 



1.1/ ^0 



|k|<X 

We thus obtained that 

Ve>0 3No = No{^r)eN: \\{A - DnMI^^^^^^^].^^^ < e, ViV>iVo, 
ending the proof. 



□ 



Remark 4.1. Let {uj\f,p]y) = {D]\f{w^), D]\f{qj\f)), and define {u,p) := {Aw,Aq). Our 
proof also shows that the field {u]\[,pn) satisfies the equation 

dtUN + {Dn o G) (V • (UJV ® UTv)) - I^AUTV + VpN = {Dn o G)(f), 

(4.9) V • = 0, 

u;v(0,x) = (DAroG)(uo)(x). 

This equation is consistent with the convergence result, since D]\j o G ^ Id, and the 
proof contains the fact that (u,p) = limjv^+oo(w7v, Qtv) is at least a distributional solution 
of the Navier-Stokes Equations (1.1). We also recall that the energy equality holds (see 
Remark 3.4), for the solution (wtvjQat) of the ADM model (4.1). 

We now prove that the solution w satisfies an "energy inequality." Observe that, from 
the previous estimates, we obtained (as a consequence of the lower bound on the operator 
-Dat) that 

WAT G ( [0 , T] ; H2 ) uniformly in N 

WAT G L°°([0,r];H2) NON uniformly in iV, 

hence obtaining an estimate for w in L°°(0,T;H2) is not trivial at all since it does not 
derive directly from the various estimates collected in the tables. 

Proposition 4.1. Let be given uq € Hq and f € L^([0, T]; Hq), and let {{wN,QN)}NeN 
be a (possibly relabelled) sequence of regular weak solutions converging to a weak solution 
(w, q) of the filtered Navier-Stokes equations. Then, w satisfies the energy inequality 

^-^i|Awf + iy\\VAwf < {i,Aw), 
in the sense of distributions, that is 

-I r \\Aw{s)f<J)'{s)ds + i^ r \\VAw{s)f<J){s)ds< r{{{s),Aw{s))(l>{s)ds, 
^ Jo Jo Jo 

for all (j) € C^(0, T) such that 4> >0 (From this one can derive the more familiar integral 
formulation (4.12). See [5, 9, 24])- This implies that w is the average of a weak (in the 
sense of Leray-Hopf) or dissipative solution u of the Navier-Stokes equation (1.1). In fact, 
the energy inequality can also be read also as 

i^lluf + z.||VAuf < (f,u). 

If we assume less regularity on the external force, as for instance f G L^([0, T]; H_i), the 
proof remains the same and we obtain the corresponding inequality 

- — IIAwlP + iy\\VAwf < (i,Aw). 
2dt - ^ ' ' 

Proof. We proved that a (relabelled) sub-sequence {wAr}„gN converges to a vector field w, 
which is solution in the sense of distributions of the filtered Navier-Stokes equations. The 
final step is to prove that the solution w satisfies the "energy inequality." In fact observe 
that the a priori estimate we proved show that 

wGL2([0,r];H2) ^ uGL2([0,r];Ho), 



and consequently the solution u belongs to L^(T3), for a.e. t G [0, T]. We need now to 
show that the L^-norm is bounded uniformly and that an energy balance holds. To this 
end we recall that the following energy equality (see Remark 3.4), is satisfied, for each 
iV > 0, by the solutions of (4.1) 

(The previous inequality has to be intended in the sense of distributions.) What is needed 
is to pass to the limit as ^ +oo. First, let us write the equality in integral form, since 
it will be most useful for our purposes: for all t £ [0, T] 

hAy'D]/\^^){t)f + u f\\VA'/'D]/\wN){s)fds 
(4.10) , 

= l\\Ay^DfMiO)f + l {A-^"Df{i),A'I^Df{^^))ds. 

Observe now that, by the estimates in (3.16-b) we have established the following conver- 
gence 

A"^d]1\^n) ^ ^1/2^1/2 (w) weakly in L2([0, T]; Hi), 
A^I'^d]1^{^n) ^^/^^^/^(w) weakly* in L°°([0, T]; Hq). 

This follows since, we are working now on a sequence {wArjTveN which is convergent in all 

the spaces previously used. In particular, we have that {A-'^/^Dj/^(wAr)}7V6N converges to 

"something," weakly in L2([0, T]; Hi), and also {A^/'^D]!j'^{^NN)}N&i has a L°°([0, T]; Hq) 
weak* limit. Finally, since we previously showed that 

Dn{^n) A-w strongly in i.^([0,r] x T3), M p < 10/3, 
D]\fiwN) ^ Aw weakly in L^([0, T]; Ho), 

the same arguments as before show also that 

^ A^/\w) weakly in L2([0, T]; Hi), 

and the uniqueness of the weak limit proves (4.11). 
Next, due to the assumptions on f we have 

^-i/2^i/2f ^ f Strongly in L2([0, T]; Hq). 

In addition, since for all N G N, WAr(O) = w(0) = u(0) € H2, and since the integral 
involving the external force f consists one term weakly converging, times another one 
strongly converging, it holds that 



lim 



lpV2<(w^)(0)f + ^\A-V2<(f),^V2^V2(^^))^^- 

= ipw(0)f + l\{,Aw)ds. 



The previous limit implies that the left-hand side of (4.10) is bounded uniformly in G N 
since 

lim sup 



< liminf 

N-t+oo 

= lim 



'\\\Ay'D]l\^^)m\' + j\A-'''Df{ilA'I^D]i\^^))ds 
\\\A^{Q)f + j\i{s\A^{s))ds. 



Next, we use the elementary inequality for the real valued sequences {aAr}„gN and {6Ar}n 



eN 



with 



lim sup Cat + liminf < limsup(aiv + &Ar), 

N-*+oo N-^+oo N^+oo 



aN:=]-\\A^I^D'-J,\^N){t)f and hN = v f \W A^l^ D'-J,\^M){s)f ds. 



(The inequality holds since we know in advance that the right-hand side is finite.) We 
infer that 

limsup-pi/2Z)]v/^(w^)(t)f + liminf / \\V A^/"^ d]1'^{v^n){s)\\^ ds 

<^\\AwiO)f + j\iis),A^{s))ds. 
By lower semi-continuity of the norm this implies that 

[ ||VAw(s)fds < liminf / \\VA^/^Dl^^WN{s)fds. 

Jo N-*+ooJq 

On the other hand, since D'^/^WN A^/'^w weakly* L~([0, T]; Hq). Again by identifica- 
tion of the weak limit we get 

\\Aw{t)f < limsupP^/2_c,i/2^^^^)||2^ 

By collecting all the estimates, we have finally proved that, for all t G [0,T], 
(4.12) ^\\Aw{t)f + ul\\VAMs)fds<^-\\Ay.{0)f + l\i{^^^ 

This can be read as the standard energy inequality for u = Aw. 

^||u(t)f ||Vu(.)f ds < ^||u(0)f + J\ns), n{s)) ds, Vt G [0, T]. 

□ 

Remark 4.2. With slightly more effort (the techniques are the usual ones as for the 
Navier-Stokes equations, see e.g., [5]) we can prove the same inequality, between to cind t, 
where to is any Lebesgue point of the function \\A'w{t)\\: For a.e. to G [0,T] it holds 

h\Aw{t)f + i^ f \\VAwis)f ds <hAw{to)f + [\tis),Aw{s))ds, VtG[to,T]. 
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